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| We investigate the subtle effects of diffuse charge on interfacial kinetics by solving the governing 

f"*) . equations for ion transport (Nernst-Planck) with realistic boundary conditions representing reaction 

kinetics (Butler- Volmer) and compact-layer capacitance (Stern) in the asymptotic limit e — Xd / L —* 
0, where \jj is the Debye screening length and L is the distance between the working and counter 
electrodes. Using the methods of singular perturbation theory, we derive the leading-order steady- 
state response to a nonzero applied current in the case of the oxidation of a neutral species into 
cations, without any supporting electrolyte. In certain parameter regimes, the theory predicts a 
reaction-limited current smaller than the classical diffusion-limited current. We also analyze the 
impedance of the electrochemical cell when a small AC current modulation is added to an applied 
DC current. At sufficiently high AC frequencies, the Maxwell displacement current is found to 
exceed the Faradaic conduction current, and experimentally observed "negative impedances" (out 
of phase AC voltage responses) are predicted close to the reaction-limited current. Overall, we 
demonstrate that the dynamics of diffuse charge plays a fundamental role in nonequilibrium surface 
reactions when the transport of one of the reacting species is coupled to the total interfacial reponse 
of the compact and diffuse layers. 
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£ , I. INTRODUCTION 

o 

The influence of diffuse-layer structure on electrochemical surface reactions has been studied for many years since 
the pioneering work of Frumkin The double layer is typically modeled as a stationary outer "diffuse layer" (or 
"Debye layer" ) of charged species in transport equilibrium and a inner "compact layer" (or "Stern layer" ) where the 
electroactive species can react at an effective electrostatic potential <& r (reaction plane) which is different from the 
electrode potential <& e by an amount A$s = $ e — <I> r , the Stern-layer voltage. In such models the only effect of the 
diffuse layer on interfacial kinetics is to modify the activation energy and effectively rescale the interfacial concentration 
of the reacting species |^-|l5|]. Within this theoretical approach, two types of models have been used: Either (i) the 
C^) • diffuse layer is treated as a continuum retaining the essential features of the mean concentration and electric field 
profiles [p] p~0[| or (ii) the interfacial double layer is treated as an idealized set of nested layers (e.g. the inner and 
outer Helmholtz planes) consisting of discrete molecular and/or ionic species [p|-^6[. When diffuse-charge effects are 
incorporated into either type of model, the diffuse layer is assumed to have a stationary, equilibrium distribution of 
charges. The dynamical effect of the diffuse layer on charge transfer kinetics has only rarely been considered since it 
was first mentioned by Levich fourty years ago 0,|l8| . 

The present study is motivated by the idea that the dynamical charging of the diffuse layer could also have an 
important effect on surface reaction kinetics in weak electrolytes (whose Debye screening length significantly exceeds 
the effective width of the compact layer) . In this work, we adopt the classical Nernst-Planck "mean- field" continuum 
model for ion transport, without considering the possible steric effects due to finite ion sizes |jg|| . We do, however, 
allow for the compact Stern layer to have a different permittivity eg from the bulk solvent ej, as described in Ref. p0[ , 
and we also consider nonlinear reaction kinetics given by the Butler- Volmer equation without any approximations 
(e.g. linearization). 

Unlike previous studies, we focus here on dynamical, non-equilibrium effects of the diffuse-layer charging on the ef- 
fective kinetic parameters measured by impedance spectroscopy. There have been few attempts in the literature aimed 
at characterizing the dynamics of the diffuse layer and its interplay with surface reactions, and existing studies have 
been mostly numerical in various approximations |^0| , ^7| , pT|j22| . Here we illustrate the power of asymptotic analysis to 
provide analytical predictions of the statics and dynamics of realistic interfacial double layers in response to both DC 
and AC currents. Asymptotic analysis has rarely been applied to electrochemical systems, and previous studies have 
only focused on static cases with oversimplified boundary conditions |^,|7|,p3[-p6| . Recently, more realistic boundary 
conditions have been proposed to model the interplay between the diffuse and compact layers, and asymptotic analysis 
has provided insight into the static response of a binary cell near the classical diffusion-limited and reaction-limited 
currents ]2Q] . In this work, the same realistic equations and boundary conditions are used to model the AC+DC 
response of a three-species electrochemical cell (composed of a binary electrolyte plus a neutral solute molecule). To 
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our knowledge, this is the first asymptotic analysis of a time-dependent, nonequilibrium electrochemical system. 

The article is organized as follows. In section [ijj, we describe the electrochemical system and present the model 
equations for ion transport and surface reactions. In section HI leading-order asymptotic approximations of steady 
solutions to these equations are derived in the case of a DC applied current. Finally in section |y] the asymptotic 
analysis is generalized to a time-dependent case by deriving the linear response to a small AC modulation added to 
a nonzero DC current. The full details of the time-dependent analysis will be presented elsewhere, but here we focus 
on the prediction of the negative impedance which has been observed in experiments and invoked as interpretation 
of the occurrence of oscillatory instabilities in electrochemical systems 



II. DESCRIPTION OF THE MODEL EQUATIONS 



A. Transport Equations 

For simplicity we consider an electrochemical cell with a planar geometry, possessing one electrode at X = 
(working electrode) and another at X — L (counter electrode), and hence we ignore the possible effect of convection 
on ion transport. (Note that we use upper-case letters to denote quantitites with dimensions to distinguish them 
from the corresponding dimensionless quantities introduced below, which are denoted by lower-case letters.) In order 
to limit our study to a small number of species, we analyze the case of a cation C which reacts at the electrodes to 
produce a neutral molecule N according to the electrochemical reaction 

C + ze-^N (1) 

in the presence of an anion A, which aids in charge neutralization. For analytical convenience, we also assume that 
the electrolyte is symmetric, zc — — za = z. 

The continuum fields representing the state of the solution are the cation, anion and neutral molecule concentrations 
(moles/volume) Cc(X,t), Ca{X,£) and Ci>f(X,t), respectively, and the electrostatic potential &(X, t). (Note that 
we use a lower case t for dimensional time, to distinguish it from the absolute temperature T.) Although we allow 
for time-dependent electric fields and currents, it can be shown that magnetic effects are negligible at the time scales 
under consideration so the electric field is given by E = —dx<&- In the dilute solution approximation, the average 
concentrations evolve according to the Nernst-Planck equations describing diffusion and electromigration 

d t Cc = Dd 2 x Cc + ^dx(C c dx^) (2a) 
d t C A = Dd 2 x C A - ^d x (C A dx^) (2b) 
d f C N = Dd 2 x C N (2c) 

where the diffusion coefficients and electrostatic mobilities of all species are assumed to be constant and equal to D 
and fx, respectively. The electrostatic potential is given by Poisson's equation 

- e 6 S* * = zF(C c - C A ) (3) 

where e& is the permittivity of the bulk solvent (assumed to be constant outside the compact layer) and F is the 
Faraday constant. 

The Faradaic current density due to the electromigration (conduction) of charged species is given by 

J F = -zF [Dd x (C c - C A ) + v(C c + C A )dx®] ■ (4) 

Although we consider time-varying electric fields which are too small to generate significant magnetic fields, we will 
show that they can be large enough to produce a non-negligible Maxwell displacement current. Therefore, the total 
current density is given by 

J(t) = J F (X, t) - e b d t dxHX, t) (5) 

which is the sum of the conduction and displacement current densities. Note that J is uniform across the cell since 
dxJ = 0, and therefore J is the experimentally controlled (or measured) current in a time-dependent situation, not 
J F . 
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B. Boundary Conditions 



At each electrode, X = and X — L, there are three boundary conditions on the ionic fluxes expressing mass 
conservation in the electrochemical surface reaction (Q) 

DdxCc + nC c d x 3> = -Jf/zF (6a) 
Dd x C A - lJ-C A d x ^> = (6b) 
Dd x C N = J F /zF (6c) 

(We adopt the sign convention that a cathodic current is negative.) In the mean-field approximation, the Faradaic 
current Jp at each electrode (which contributes to the surface reaction, unlike the displacement current) is related to 
the local concentrations and potential through the Butler- Volmer kinetic equation (applied at X — and X = L) 

J F = -K C C exp (-a zFA<$> s /RT) + K R C N exp {a R zFA$ s /RT) (7) 

where Ko and K R are the oxidation and reduction kinetic rate contants, respectively, Cc and Cm are the interfacial 
concentrations, ao and a R are the transfer coefficients (which are set equal to 1/2 below) and A<f>s = $ e — $ r is the 
voltage drop across Stern's compact layer mentioned above. Following Frumkin we imagine applying (R) at the 
outer Hclmholtz plane, or any other convenient molecular distance which acts as the edge of the continuum region. 
Note we do not apply the Butler- Volmer equation across the entire interface (including both the compact and diffuse 
layers) as is frequently done. 

No other boundary conditions are typically mentioned in electrochemistry textbooks because none are needed when 
the common assumption of electroneutrality is made |3l| . In this work, however, since we treat diffuse charge explicitly 
we need another boundary condition on the potential. Ignorance of this boundary condition is typically hidden in the 
"zeta-potential" £ (the potential drop across the diffuse layer A&p,) or equivalently the total charge in the diffuse layer. 
In colloidal science these quantities can be assumed to be constant properties of a surface, but in electrochemistry, the 
zeta potential of an electrode should be determined self-consistently from the microscopic electrochemical boundary 
conditions. We will see that the zeta potential of a working electrode can vary widely with electrochemical conditions, 
especially in time-dependent situations. 

A general expression for the missing boundary condition has recently been proposed based on the Grahame model of 
the interface, which contains the Stern model as a special case ^TJ. Here we adopt the Stern model, which postulates 
a constant capacitance Cs for the compact layer 

A$5 = 1 esdx^/Cs at X = L ^ 

where es is an effective permittivity for the Stern layer. Note that we do not consider, in this model, any chemical 
adsorption effect. In this work, we model experiments in which the current J(t) is prescribed and the voltage difference 
between the working electrode and the reference electrode A<i> tot = <i> e — 3> re / is determined by solving the equations. 
Note that 

As - es/C s (9) 

defines an effective width for the Stern layer. We do not claim that As corresponds to any well-defined molecular 
distance, but only that it captures the combined effect of the compact layer capacitance and permittivity. In our 
numerical calculations below, we assume that the permittivity of the Stern layer is es — 10e o , compared with the value 
£b = 80e o for a water solvent, consistent with previous studies f32fl . From the equation As = es/Cs and the literature 
data on Cs ~ 80/^F.cm -2 |32| an estimate of As ~ lA can be computed. We also assume that As is independent of 
the local ionic concentrations. 



C. Dimensionless Equations and Boundary Conditions 



The first step in any asymptotic analysis is to scale all quantities appropriately and identify the relevant dimen- 
sionless groups. Here we scale length to the electrode separation L, time to the diffusion time across the cell L 2 /D, 
potential to the thermal voltage zFD/RT, concentrations to the mean anion concentration C* 



lJo 



C* = - I C A {X)dX (10) 
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and current to the classical diffusion-limited current of Nernst 

AzFDC* , s 

Jdl = ^ . (11) 

The dimensionless variables at these scales are 

x = X/L, T = tD/L 2 , j = J/J L , Ci = Ci/C* (i = C,A,N), <p = zF<S>/RT. (12) 
We also introduce the mean concentration of charged species c(x,t) and the charge density p(x,t) 

c = i(c c + c A ) (13a) 

P = 7}(cc ~ c A ) (13b) 

In terms of these variables, the Nernst-Planck and Poisson equations take the dimensionless forms (using the Einstein 
relation fjt = zFD/RT) 

d T c = d 2 x c + d x {pd x dp) (14a) 
d T p = dlp+d x {cd x <l>) (14b) 
9 t cn — d 2 CN (14c) 
-e 2 dl<}> = p (14d) 

where the dimensionless group e = \d / L indicates that the natural scale for charge screening is the Debye length \d 
given by 

2 e b RT 

= 2F^ (15) 

rather than the electrode separation L. In this work we consider the typical case in which e is much smaller than 
unity (for a cell of typical length L = lcm and < lOOnm, e < 10~ 5 ), which is the basis for our asymptotic analysis. 
Since e multiplies the highest derivative in the equations, it is a singular perturbation which must be treated using 
the methods of boundary layer theory |}3| . Note that the expression for the dimensionless current (which is uniform 
across the cell since dj/dx = 0) 

j(t)=j F (x,t)+j d (x,t) (16) 

where 

3f = -~{d x p + cd x (f>) (17a) 

3d = -^ 2 d T d x 0(x,t) (17b) 

involves the perturbation parameter multiplying the displacement current density jd, which indicates that it only 
becomes important at high frequencies (and/or small length scales). 

The flux boundary conditions at x = and x = 1 take the dimensionless forms 

d x c + pd x (p = —2jp (18a) 
d x p + cd x <p = -2j F (18b) 
d x c N = ij F (18c) 

where the dimensionless Faradaic current at each electrode is given by 

jp = -ko (c + p) exp (~a o A0 s ) + k R c N exp (a R A(f> s ) (19) 

where ko = KoL/AD and k R = KrL/AD are the dimensionless kinetic constants and Atfis — zFA^s/RT. The 
dimensionless Stern boundary conditions 
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involve a second dimensionless group 6 — Xs/^d, which is the ratio of the effective widths of the compact and diffuse 
parts of the interfacial double layer. The Gouy-Chapman model (with no compact layer) corresponds to the limit 
(5 = while the Helmholtz model (with no diffuse layer) corresponds to the limit 6 — > oo. In this work, we consider 
the full range of values for 8, but always require K<1. 

It is ubiquitous in electrochemistry to assume that the total charge in the cell 

q(t) = f p(x,t)dx (21) 
Jo 

is identically zero, but this "macroscopic intuition" (which is valid in the limit e — > 0) is flawed whenever diffuse 
charge is being treated explicitly (e > 0), as we do here, because there is a total charge (of order e) in the diffuse 
layers which exists to satisfy the nonlinear and asymmetric reaction conditions at the electrodes |2C| ] . 

The fact that q(t) can vary in time does not contradict the conservation of law for charge because dq/dt is related 
by (14b) to the Faradaic currents at the electrodes 

d T q = f d T p(x,t)dx = 2\j F (0,T)-j F {l,T)} (22) 
Jo 

rather than the total currents (which are equal). Since the total current j = jp + jd is uniform across the cell, any 
changes in the total charge of the cell are caused by temporary differences between the displacement currents at the 
two electrodes. 



III. STEADY STATE RESPONSE TO A DC CURRENT 



A. The Smyrl-Newman Equation 



In this section, we derive leading-order steady-state solutions of ( ff4| ) and (p^)~(p0|) in response to a constant applied 
DC current density j , which are valid after transients have died away, following the analysis of Ref. pp| . Integrating 
the steady Nernst-Planck equations once, we arrive at a system of ordinary differential equations 



d x c + pd x (j) = -2j 
d x p + cd x (f> = -2j 
d x c N = Aj 
-e 2 2 x( f> = p 

By subtracting the first two equations we obtain 

d x (c - p) = {c- p)d x 4> 



(23a) 
(23b) 
(23c) 
(23d) 



(24) 



which is easily integrated to show that (in steady-state) the anions are in Boltzmann thermal equilibrium throughout 
the cell 



ca = c — p cx exp(</>) = exp(— zaF<& / RT) 



(25) 



since they do not react at the electrodes. 

Our analysis begins by combining the three coupled equations (23a), (23b) and (23d) into a single equation for the 
potential, which is decoupled from the fourt h equ ation (23c) for the neutral species concentration (except through 
the boundary conditions). Substituting (23d) into (23a) and integrating we obtain 



e 2 
c(x) =c - 2jx + — (d x (f)) 



(26) 



where c G is an integration constant to be determined by the boundary conditions. If we substitute this expression 
and (23d) into (23b) we obtain a third order nonlinear equation for the potential (2(J 



+ (2jx - c ) d x (j) = 2j 



(27) 
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A similar equation was first derived by Smyrl and Newman for the electric field of a rotating disk electrode near the 
diffusion-limited current (7). O nce t his equation is solved for the potential <j>, the concentration c and charge density 
p are computed from (|2^) and fl23d|) , respectively. 

The Smyrl-Newman equation (|27|) cannot be solved analytically in terms of elementary functions, but due to the 
singular perturbations (terms involving e), it is even difficult to solve numerically. Both of these difficulties can be 
conveniently resolved with asymptotic analysis [fU?]!^] . The idea is to construct a uniformly valid approximation for 
all x (in the limit e — > 0) by adding the leading order solution to the "outer problem" in the neutral bulk region 
e < a; <C 1 - e to the leading order solutions to the "inner problems" in the charged diffuse layers < x <C e and 
1 — e <C x < 1 and subtracting the overlapping parts through asymptotic matching, which is valid well below the 
diffusion- limited current \j\ = 1. As emphasized by Newman ||, the asymptotic approach used here reduces to the 
familiar Gouy-Chapman (Poisson-Boltzmann) theory of the diffuse layer in the special case 6 = 0, but it is more easily 
generalized to the time-dependent case. We also study for the first time the novel nonlinear effects of the coupled 
Stern and Butler- Volmer boundary conditions when 6 > 0, which include the possibility of a reaction-limited current. 



B. Leading Order Solution of the Outer Problem 

In the outer region, we set e = in (E7h to derive the leading order approximations 



d x ct> = 2j/(2jx-c ) + 0(e) (28a) 
c{x) = c -2jx + 0(e) (28b) 
p{x) = 0(e 2 ) (28c) 



from (J26|) and (23d), where the integration constant c D can be evaluated using the conservation of anions over the cell 



Jp 1 CA(x,t)dx — 1 at zeroth order in e. 

c = l+j (29) 

since q(t) = 0(e). Note that the outer region is neutral (p = 0) at leading order, as is commonly assumed in 
electrochemistry. The bulk potential profile at leading order is obtained by a simple integration 

<t>{x) -c/>e = + log [( Co - 2jx)/c ] + 0(e) (30) 

where A(j>i = <f> e — (f)(0) is the total (leading-order) potential drop across the interface consisting of contributions from 
the diffuse (Debye) and compact (Stern) layers 

= Acf>D + A<j> s (31) 

which must be calculated by asymptotic matching with the inner problem. (Note again that A<pr> is the familiar zeta 
potential.) The equation for the neutral species concentration is also easily integrated, writing the conservation over 
the cell of the sum of the cation and neutral species numbers at zeroth order in e. 

c N = 7 - 2j(l - 2x) + 0(e). (32) 

Typical bulk concentration, electric field and potential profiles are shown in Fig. ^. (We set 7 = 1 in all the numerical 
results presented here.) The potential profile <p is slightly curved (logarithmic dependence on x), and the electric field 
d x 4> is not constant since the migration term is a nonlinear function of the concentration c (which varies linearly in 
the bulk when j ^ 0). The profiles of cf> and d x (f> corresponding to opposite values of j (e.g. 0.4 and -0.4) are not 
symmetric across x-axis because the product of the reduction of cations is neutral. On the other hand, the profiles of 
c and cjv are symmetric when the current is reversed. 



C. Leading Order Solution of the Inner Problem 



The leading order approximation for the diffuse layer potential at the x = electrode is obtained by changing 
variables to the inner coordinate y — x/e in fl27 



- dl<t> + ~ (d y 0) 3 + c d y 4> = O(e). (33) 



G 



where we ignore the O(e) terms. Multiplying by 2 dy<j> and integrating (using the matching conditions d y 4>, dy 
as y — > oo ) we obtain 



d 2 y 4> = -dy^co + (9^)74 (34) 

which is a first-order separable equation for d y 4>. The next integration can be performed with an hyperbolic substitu- 
tion d y 4> = ±2y/c^csch(u) (where csch(it) = l/sinh(u)), which yields the trivial equation d y u = -Jc^, whose solution 
is u = ^fcZ(y + K) for some constant K, Therefore, the electric field in the diffuse layer is 

d y (j } = ±2^r o csch{^r {y + K))+0{t) (35) 

where K is related to the Stern-layer voltage through the boundary condition at y = 

Acb s = T2S^csch{^K). (36) 

The upper sign refers to potentials below the potential of zero charge ("p.z.c") with A4>d,^4>s < 0, and the lower 
sign to potentials above p.z.c. 

Transforming to the inner coordinate in ( p6| ) and (23d) 

c = c + - {d y dpf + 0(e) (37a) 

p = -8^ + 0(e) (37b) 

and substituting ( |35| ) for the electric field, we obtain the leading order concentration and charge density profiles 

c = c [l + 2csch 2 + K))} + 0(e) (38a) 

p = ±2c csch(^{y + K))coth(^/c~{y + K)) + 0(e) (38b) 

The constant of integration K, or equivalently the Stern layer voltage A0s, is determined by solving a transcendental 
equation provided by the Butler- Volmer kinetic boundary condition at y = (see below) 

j = -ko(c + p) exp (-a o A0 s ) + &r(7 - 2j) exp (a R A<p s ) (39) 



where c(0) + p(0) = cc(0) is a (generalized) Frumkin correction |31|. 

The expression for the electric field j35|) can be integrated analytically to obtain the diffuse layer potential profile 
at leading order 

e - (j) = A(f>i + 4tanlT 1 f e -^ (y+K A (40) 



which is essentially the Gouy-Chapman solution to the Poisson-Boltzmann equation, although our boundary conditions 
(specifying K) are substantially different from the classical theory, which corresponds to the limit 8 — > 0. To clarify 
this connection, we let ip — Afa + <j>, and observe (following some algebra) that the solution ([lo]) satifies 

dytp = 2V^sinh(W2) (41) 

and 

dfa = c sinh(^) (42) 

which is the Poisson-Boltzmann equation. (Changing the current simply changes the bulk concentration c = 1 + j.) 
The Smyrl- Newman equation ( |27j ) is more general than the Poisson-Boltzmann equation, but the former is equivalent 
to the latter below the diffusion- limited current J2(J. Combining ( |36| ) and (41), we obtain a relation between the Stern 
and diffuse layer voltages 

A(j) S = 2Sy/c~smh{A(j) D /2) (43) 
and from (|I|) the constant K can be related to the diffuse layer voltage 

tanh(A</> jD /4) = e -^ K . (44) 
To determine K we must solve the transcendental system (|3^) and (§91) numerically. 
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D. Application of the Electrochemical Boundary Conditions 



The behavior of K as a function of the ratio of the kinetics constants Rk = kfi/ko and the current density j is 
reported in Fig. || (with kn — 1 held constant). From ([36]) it is clear that large values of K correspond to small values 
of A(f>s when j = 0, i.e. the symmetric case c G = 1. The introduction of a nonzero value of j breaks the symmetry of 
the K versus log(i?fc) curve. The plot of K versus j in Fig. |](b) illustrates the dissymmetry of the redox couple, since 
when the current is positive the oxidation of neutral species produces positively charged cations which accumulate 
in the diffuse layer, thus hindering the negative charging of the diffuse-layer. In Fig. |^ we also display the Stern and 
total interfacial voltages, A^>s and A</>j = A0s + A<pD- We observe that when j increases beyond 0.3, A(j>s and Afa 
reach arbitrarily large values to compensate this excess of cations (see below). For j — 0, the plots of A4>s (Fig. §(c)) 
and A(j>i (Fig. |](e)) are symmetric with respect to Rk = 1. This behavior is confirmed by the plot of the different 
spatial profiles c, p, d y cj> and (p, for two different values of j, in Figs. || and || 

In Fig. H we focus on the case of zero current. When the kinetic constants are equal Rk = 1, this correspond to the 
p.z.c, but when Rk ^ 1, the system builds a charged diffuse layer to equilibrate the oxidation and reduction currents 
on each electrode. For comparison, we have also added in (dashed lines) the stationary Poisson-Boltzmann solution 
for j = 0, which corresponds to choosing K such that 5 = in ([h)|). 

In Fig. ^, a positive current j — 0.3 has been chosen, with symmetric kinetic constants (Rk = 1). In this case, 
the system builds a charged diffuse layer to favor either reduction (negative Faradaic current) or oxidation (positive 
Faradaic current). In each figure, two values of Xd have been chosen, one corresponding to a weak electrolyte 
concentration Xd = 100As, where the diffuse layer has a strong effect on the charge distribution and potential 
profiles, and another corresponding to a strong electrolyte concentration with Xd = As, where there is a much more 
significant potential drop in the Stern layer (and hence less influence from the diffuse layer). We note that for j = 0.3 
the charge and electric field profiles shown in Fig. [I] start to diverge at y = 0, particularly so at the small value of 
5 = 1/100. 

To better understand this unexpected behavior observed for j > 0.3, we plot the influence of the ratio Ad/As = S^ 1 
on the potential drops of the diffuse layer A(j>D, Stern's layer A<f>g and the whole double layer A<^ = A(f>D + A(f>s in 
Fig. ||. We notice that an exchange of potential drop between Stern's layer and the diffuse layer occurs when the size 
of the diffuse layer becomes larger than As . The potential drops in Stern's layer and diffuse layer compensate exactly 
when j — 0, whereas the potential drop in the diffuse layer diverges when Xd increases at a nonzero current larger 
than 0.3. 



E. Theoretical Polarograms 



For the electrochemist, polarogramms bear a great significance. In the present study of stationary regimes in finite 
size cells such a representation of current j versus the total voltage between the reference and working electrodes 
Atfitot = 4>e — 4>ref is shown in Fig. || This includes leading order voltage drops calculated for both double layers 
and the bulk electroneutral region. It confirms again than for j positive large than 0.3, a limitation of the current 
occurs which is very different from diffusion limited current which are usually encountered in electrochemical cells. 
The standard current-potential curve ("PB", plotted as a dotted line in Fig. ||) shows a limiting current equal to 7/2 
which occurs when cm — ► at the working electrode surface. When j > the production of cations increases due to 
oxidation of neutral species, to maintain the imposed current, the system usually increases A0s to compensate this 
excess. When S — Xs/Xp becomes much smaller than 1, the potential drop in the diffuse layer increases much more 
rapidly than the potential drop in the Stern layer, there is an enhancement of the electric field in the diffuse layer 
which tends to strip the cations from the anode. For a given current density (j > 0.3 in this case) the interfacial 
cation concentration tends to zero. This limitation of the current, however, is not due to transport but rather the 
reaction kinetics coupled to the diffuse layer charging. 

The significance of j — 0.3 in our numerical results is most easily understood in the limit of a very weak electrolyte 
with Xd S> As (5 — * 0) where it corresponds to a "reaction-limited current" , generally given by 

JRL = (45) 
1 + 2k R 

whenever this value is smaller than the Nernst diffusion-limited current 7/2, when c/v(0) = 0. (Recall that we have 
chosen fcfj ; = 7 = 1, so that j^h = 1/3 and 7/2 = 1/2.) Instead of being caused simply by transport limitations, the 
reaction-limited current jm, is determined by the electrochemical boundary conditions. Note that when 6 = (no 
Stern layer), and the Butler- Volmer equation ([39]) takes the simpler form 
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j = ~ko{c + p) + k R (j - 2j), 



(46) 



and in this case, a limitation of positive current occurs because the first term representing the cathodic current is 
negative (but small) while the positive second term representing the anodic current decreases with applied current as 
the concentration of neutral species is reduced (but not necessarily depleted) at the electode. When 5 > (however 



small), the situation is quite different because the system can avoid the reaction limitation in (39) by producing a 
large positive Stern-layer voltage A<frs 

A0 S ~ — log ( 3 ) (47) 
an. \k R {j-2j)J 

which suppresses the negative cathodic current and exponentially enhances the positive anodic current enough to 
maintain the applied total current. This also produces a large diffuse layer voltage for j > jrl according to fl4S|), 
although this effect is reduced as 5 is increased. 

The actual limiting current in the system is controlled by the smaller of the reaction-limited current (^|) and 
the diffusion limited current \j\ = 7/2. Note that for negative currents one observes a diffusion limiting current at 
j = —7/2 which corresponds to the annulation of the neutral species concentration at the counter electrode. Therefore, 
the system never reaches the classical limiting value j = — 1 which would be obtained with a strongly supporting 
electrolyte (purely diffusing system). 

IV. ANALYSIS OF THE LINEAR RESPONSE TO AN AC CURRENT MODULATION 

In this section we derive the leading-order response of our simple electrochemical system to small amplitude AC 
current perturbations (around a steady nonzero DC current) using the same mathematical methods as in previous 
section, now applied to the Laplace transforms of the linearized equations. There have been very few attempts in 
the existing literature |34|-|36|] to compute the impedance spectra for electrochemical systems with weak electrolytes 
(and large Debye lengths) in the absence of a supporting electrolyte. The modeling of voltammetric curves has 
received more attention p7|-|4l[] , although most of these studies either neglect the diffuse layer entirely or assume its 
ionic concentration profiles are those of equilibrium. In this section we show that the diffuse-layer cannot always be 
approximated by an additional fixed capacitance in parallel with the Stern-layer capacitance because its dynamics 
can strongly influence the net interfacial impedance at intermediate frequencies. 

To compute the response of the electrochemical cell to small applied AC current modulation 5j of a DC current j 

j=3 + Sj (48) 
we add a corresponding small modulation to each of the variables 



c 


= c + 


5c 


(49a) 


Cjv 


= Cjv 


+ Scn 


(49b) 


p 


= p + 


5 P 


(49c) 




= «H 


-S(f> 


(49d) 








(49e) 



about the leading-order approximation of the steady-state solution described above (denoted by an overline), and 
linearize the time-dependent model of equations (|l4j) around the stationary state to obtain the following system: 

d T Sc = d 2 Jc + d x (pd x S4> + Spd x 4) (50a) 

d T 8p = dl8p + d x (8cd x 4> + cd x 5(j)) (50b) 

d T 5c N = d 2 x 5c N (50c) 

5p = -e 2 d 2 J4> (50d) 

The electrochemical boundary conditions at both electrodes (x = 0, 1) are also linearized about the steady-state 

— 2SjF = d x 5c + pd x 5(f> + 8pd x 4> (51a) 
—25 jF = d x 5p + 5cd x (j) + cd x 5<p (51b) 
A5j F = d x 5c N (51c) 
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SjF = [fcoao(c + p)exp(-a A</)s) + k R a R c N exp(a R A^ s )] 5A(j> s (51d) 
-k (Sc + 5p) exp(-a o A0 s ) + k R 5c N exp(a R A<j>s) (51e) 



where the Stern-layer voltage is decomposed as the other variables A(f>s = Acf>s + SA4>s- In the time-dependent 
setting, the current modulation is broken into separate contributions from the Faradaic conduction and displacement 
current densities 

5j = S JF + Sj d (52) 

where 

Sj F = -^(d x 5p + 5cd x <j> + cd x 8<j)) (53a) 

Sj d = ~e 2 d T d x 5<t> (53b) 

The unperturbed DC current is entirely Faradaic (j = jp, jd = 0), but the displacement current becomes important 
in the AC modulation, especially at large frequencies and small length scales. 

To compute the impedance of the electrochemical cell, we consider a small amplitude AC perturbations of the form 

Sj = jo sin(wf) (54) 

and take the Laplace transform in time of the system of transport equations (|5C| ) 

sSc = d x Sc + d x (pd x S~4> + Spd x (f>j (55a) 
sSp = dl5p + d x Ucd x & + cd x 6~<i?) (55b) 

sSc n = d x Sc n (55c) 
Sp = -e 2 d 2 J4> (55d) 

(55e) 

where s = iujL 2 / D is the Laplace variable conjugate to time (the dimensionless forcing frequency, or imaginary decay 
rate) and tilde accents are used denote the transformed variables 6f = CSf. Since they do not depend explicitly 
on time, the interfacial boundary conditions ( pT| ) are the same for the transformed variables as for the original ones. 
Note that the displacement current modulation takes the form 

Sjd = -\e 2 sd x 5~4> (56) 

in Laplace space. Once we solve these equations for the perturbed variables, the (complex) impedance of the cell is 
calculated as the ratio of the voltage response versus the current modulation in Laplace space 

5A~(j) tot 

where 8A4> t ot is the total voltage modulation (in Laplace space) between the working electrode and the reference 
electrode. 

The analysis begins by combining the transport equations (^5|) following the same steps as in the derivation of the 
Smyrl-Newman equation above, which leads to a system of two linear equations for 5<f> and 8c: 

-se 2 d 2 x 8~(t> = -e 2 d*S~<p + d x Ucd x $+cd x 8~<t?) (58a) 
sS~c = d 2 Jc- e 2 d 2 x [d x 4:d x 5(j^ (58b) 

Although these equations are linear, they cannot be solved analytically in terms of elementary functions, and they 
also possess singular perturbations which hinder numerical progress. As in the steady-state case, asymptotic analysis 
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in the limit e — > once again allows us to solve the leading-order problem by decomposing the cell into a bulk 
electroneutral region which must be asymptotically matched with charged diffuse boundary layers. The resulting 
leading order equations for the inner and outer problems are difficult to solve analytically in this case, but they present 
no trouble for a numerical solution because the singular perturbations have been removed analytically. The details 
of these calculations will be presented elsewhere, but here we present a simple example of the rich time-dependent 
electrochemical behavior that can be predicted by the mathematical model. 

In the framework of the asymptotic analysis, the total impedance of the cell is conveniently broken into separate 
leading-order contributions from the various regions 

Z{ U ) = ^k+j^g+j^ (59 ) 

where SAcfrs, SA(j)o and 5A(j)b are the voltage modulations in the Stern layer, the diffuse layer and the bulk electroneu- 
tral zone, respectively. The asymptotic analysis for the limit e — > assumes that the diffuse layers are stationary 
(quasi-steady) at leading order, which is valid when the forcing frequency is below the Debye frequency 



D 



or in terms of the dimensionless frequency 



uL 2 1 . , 

This does not limit the applicability of our analysis because impedance spectroscopy experiments are conducted well 
below the Debye frequency. 

It is interesting to consider the displacement current density over this range of time scales. From ([36]) we see that 
when \s\ <C 1/e 2 or e 2 s <C 1 the displacement current Sjd is always negligible in the outer electroneutral region. 
However, the situation is very different in the diffuse layers, where the natural length scale Xd is much smaller. In 
the diffuse layer, we write the displacement current in terms of the inner coordinate y — x/e, 



S Jd = --sedydcj) (62) 



and observe that it becomes comparable to the Faradaic current when |s|A_d<5A0s/As — 0(e 1 ). Since the kinetic 
constants (ko, &r) are included in the value of 5A05, they will have strong influence on the current displacement and 
on the "crossover" frequency where the displacement and Faradaic currents become comparable. 

Note that for the case of non-reacting system (zero kinetic constants, and zero Faradaic current) the characteristic 
frequency of Debye layer charging has recently been derived in the completely different context of electrohydrodynamic 
response to a spatially periodic applied AC voltage Q : 

" C= A^L ■ (63) 

This prediction is recovered by our model in the limit where Xjj 3> A5. 

In Fig |?], a plot of the two current modulations versus the AC frequency u gives a clear insight into the interchange 
of the reaction and diffuse layer charging processes with frequency. (The amplitude of the stationary current in 
this example is j = 0.3, near the steady-state reaction-limited current Jbl)- The crossover frequency where the 
displacement and Faradaic current amplitudes are comparable is strongly affected by the parameter <5 = Xs/Xd- 
The larger is 5, the smaller is the diffuse layer size and the higher the frequency at which the displacement currents 
dominate, given (very roughly) by (T 1 . 

The arguments of the current modulations are such that when Xd increases the phases become shifted by tt, showing 
an out of phase AC voltage response. This phenomena is often termed "negative impedance" p7|-^o| and has already 
been reported in the electrochemical literature as the source of instabilities (stationary or oscillatory). The impedance 
spectra corresponding to the current modulations reported in Fig. Q are shown in Fig. |^. The two examples 5=1, 
e = 1CP 8 and S — 0.01, e — 10 -6 are consistent with the phase shift that we have described in the current modulation 
curves. The amplitude of the impedance also increases very dramatically when the current density goes beyond 0.33. 
As already discussed in the stationary case, when the system reaches this "reaction-limited current" situation, a 
strong potential drop is established in the diffuse layer, which can explain the amplification of the electrochemical cell 
impedance. We hope to elaborate on this aspect of the theoretical predictions in comparison with real experiments 
in a future communication. 
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V. CONCLUSION 



We have performed a leading-order asymptotic analysis of a realistic set of transport equations and boundary 
conditions for a three-species electrochemical cell in which cations are reduced into a neutral molecules in the presence 
of anions, in both stationary and time-dependent nonequilibrium situations. By solving the full nonlinear problem 
at leading order (numerically where necessary) , we have made the first systematic exploration of the influence of the 
parameter 5 — \s/\d (introduced in Ref. |2fJ] to describe the coupling of the Stern and diffuse layers) on stationary 
and time-dependent interfacial kinetics. The mathematical model is able to predict very subtle effects related to the 
coupling between diffuse charge and interfacial reactions, such as the existence of both reaction-limited and diffusion- 
limited currents and negative impedances with nonnegligible displacement currents are large frequencies. The strength 
and novelty of these theoretical predictions is their firm analytical basis starting from realistic transport equations 
and boundary conditions. 

An interesting extension of this study would be to electrochemical reactions involving the reduction of anions (or 
the oxidation of cations). Such systems have already been the subject of many investigations, both theoretical and 
experimental, but most these studies have not considered simultanously the dynamics of the diffusion and diffuse 
layers and their possible interplay during the interfacial reactions, as we have in this work. Such situations would 
complicate the present analysis signficantly because there would be more than two charged species, and it appears 
that any analytical solution of the inner problem at leading order would not be possible in such cases. However, as 
long as e is small (Ad <C L), the mathematical machinery of asymptotic analysis can be used to derive simplified, 
well-behaved equations at leading order, which would be straightforward to at least integrate numerically. 

Even restricting ourselves to the simple electrochemical system studied here, there are number of interesting direc- 
tions for future work. Higher order corrections in e could be derived (with some numerical computation required). It 
would also be interesting to solve the full system of equations with e = O(l), where the asymptotic analysis breaks 
down. This limit, which corresponds to very small cells on the order of the Debye length, has increasing importance 
as microelectrochemical systems reach smaller and smaller length scales. It would also be interesting to solve the AC 
current problem more globally in the nonlinear regime, where our local linear stability analysis is no longer valid. 
In order to study nonlinear responses, such as those encountered in voltammetry or chronoamperometric methods, 
the full nonlinear system of equations would have to be solved numerically. Nevertheless, if e is small (as is usually 
the case) asymptotic analysis provides the appropriate starting point for the mathematical modeling of diffuse-charge 
effects in electrochemical systems. 
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VI. APPENDIX A. LIST OF SYMBOLS 



D diffusion coefficient 

H mobility 

zc,za,z charge numbers zc — — za = z 

F Faraday's constant 

R universal gas constant 

T absolute temperature 

t time 

Ko,Kr kinetic constants 

ko,kn dimensionlcss kinetic constants 

Cs capacitance of the Stern layer 

6b permittivity of the bulk solvent 

£5 effective permittivity of the Stern layer 

As effective width of the Stern layer 

\d Debye screening length 

L distance between the electrodes (working and counter) 

S =^s/^D 

e =X D /L 

J current density 

Jf Faradaic current 

Jdl diffusion-limited current 

Jul reaction-limited current 

Cc, Ca, Cn concentrations of the species C, A and N 

$ electric scalar potential 

C* mean anion concentration 

cc,ca,cn dimensionlcss concentrations 

c dimensionlcss average concentration of charged species 

p dimensionlcss charge density 

<j) dimcnsionless electrical potential 

j dimensionless current density 

jf dimensionlcss Faradaic current 

jd dimensionlcss displacement current 

A<fii dimensionless double layer voltage 

A(f>D dimensionless diffuse layer voltage 

Ac/)s dimensionless Stern layer voltage 
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FIG. 1. Steady-state profiles of <f> (a), dtfr/dy (b), c (c), and cjv (d) at leading order in the "outer" (bulk electroneutral) 
region for j = 0.4, and —0.4 and 7 = 1. ao = Qfl = 1/2 
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FIG. 2. (a) K versus R k = k r /k in lin-log representation, (j — 0). (b) K versus j for = 1. (c) A0 S versus Rk in lin-log 
representation, (j = 0). (d) A^>s versus j for R k = 1. (c) Acfn versus R k in lin-log representation, (j = 0), 7 = 1. (d) A(f>i 
versus j for R k — 1. Two different ratios 8 are illustrated in the six panels: circles S = 0.01, triangles 5=1. ao = ctR = 1/2 
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FIG. 3. Steady-state profiles of c (a), p (b), d y (p (c), <j> (d), cc (e) and ca (f) at leading order in the "inner" (diffuse) layer 
at y = for j = 0, R k = 5, (fc r = 5 and k D = 1), K = 0.972 for S = 0.01 (circles) and K = 1.629 for 5 = 1 (triangles), For 
comparison the Poisson-Boltzmann's profiles are reported as dashed lines which correspond to the limit 5 = 0. 
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FIG. 4. Steady-state profiles of c (a), p (b), d v (j> (c), cj> (d), cc (e) and (f) in the diffuse layer at leading order for j = 0.3, 
R k = 1, K = 0.523 for 8 = 0.01 (circles) and K = 1.394 for 6 = 1 (triangles). 
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FIG. 5. Dependence of the potential differences A(f>D (diffuse layer), A0s (Stern's layer), A<f>i = A^u + A(j>s (double layer) 
with — log 5. The different parameters Rk and j corresponding to each plot are inserted in panel (a). 
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FIG. 6. Polarograms obtained from the asymptotic analysis for (a) kn = ko — 1 and (b) ka — ko = 0.01. 
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FIG. 7. The displacement and Faradaic current modulations versus the dimensionless frequency coL 2 /D. Results for two 
different ratios 8 — A^/Ad are shown, (a) log of the modulus of current modulation Sj. (b) argument of the current modulation 
8j. As = lA and L = 1 cm, so Ad = lOnm corresponds to 5 = 0.01 and e = 10~ 6 , while Ad = lA corresponds to 8 = 1 and 
e = 10~ 8 . Other parameters are D = 10 _5 cm 2 s _1 , j = 0.3, fcjj = ko = 1, 7 = 1, ao = &r = 1/2. 
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FIG. 8. Impedance spectra obtained from a linear stability analysis around the steady state j = 0.3. (a) Nyquist plot, (b) 
and (c) Bode plots for the modulus and argument of the total impedance Z(s — iuo) = A(f>tot(s)/Sj(s). Same parameters as in 

Fig. m 
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